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The inclusive electromagnetic responses in the quasi-elastic region are calculated with a model 
which considers the terms of the cluster expansion containing a single correlation line. The validity 
Qs ' of this model is studied by comparing, in nuclear matter, its results with those of a complete 

O^ . calculation. Results in finite nuclei for both one- and two-nucleon emission are presented. 

^-> ■ I. INTRODUCTION 

CJ ■ 

^^ ' Aim of the nuclear many-body theories is the calculation of the properties of finite nuclear systems by having as 

00 , only input the bare nucleon-nucleon interaction. The approach we follow to solve the nuclear many-body problem is 

based upon the variational principle 

T-H ■ 

o 
o 



The solution of eq. (y) is exactly that of the Schrodinger equation if there are no limitations on the Hilbert space 
spanned to search for the minimum of the energy functional. In our approach we search for the minimum within a 
Q\ _ subspace composed by states of the form 
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|*o)=^|$o) , (2) 



where |$o) is a Slater determinant, formed by a set of orthonormal single particle wave functions, and _F is a many- 
O , body correlation function. 
2 The quantities to be compared to the experiment are calculated by evaluating the mean value of the corresponding 

operators between the states of eq. (0) fixed at the minimum of the energy functional: 



^: ^^^" {L\F+F\<i>o) ' ^^^ 

■ ■ ' where we have indicated with Q a generic operator associated to an observable. 

A good variational ansatz for realistic calculations requires that the correlation function F has the same operatorial 
dependence of the Hamiltonian. In this report, however, we assume that the many-body correlation is a scalar 
function which can be written as a product of two-body correlation functions: 

A 

F{r,,r2,....,rA,)=llf{nj) (4) 

i<j 

where r-y is the distance between the nucleons i and j. The choice (Q) of the correlation function leads to the following 
expression for the energy functional: 
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($o|g[l + fe(ri2)][l + hjns)] [1 + fe(r23)]...|^o) 
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In the above equation we have used the function h{rij) = /^(?'y ) — 1 to show the mechanism of the cluster expansion. 
The simplest term is the mean value of the hamiltonian between two uncorrelated states, which is obtained by 
picking up all the "ones" in the various terms in eq. (0). More complicated terms are obtained considering a single 
hijij) function multiplying the "ones". The procedure continues considering terms with two, three, etc. h(rij) 
functions. We make a classification of the various terms based upon their topology. With the Fermi Hypernetted 
Chain (FHNC) technology we calculate all the terms of a certain type. Details of the application of the FHNC to the 
description of the ground state of doubly closed shell nuclei can be found in refs. [Q— [Q. The recent progresses of 
these calculations, done with state dependent correlations, are presented in another contribution to this conference 



II. THE MODEL: GROUND STATE 



In the theoretical framework wc have presented, the evaluation of an operator mean value, eq. (|^), requires the use 
of the full FHNC computational scheme. For the evaluation of the one-body density distribution, defined as 
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we have developed a simplified model ||6|,[7| . After performing the cluster expansion of both numerator and denominator 
to eliminate the so-called unlinked diagrams, we retain the terms with a single correlation function h. The expression 
for the density obtained in this approximation is: 
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where the subindex L indicates that only the linked diagrams have to be calculated. The diagrams representing the 
terms considered in our model are shown in fig. nl In this figure the white circles represent the point ri which is not 
integrated, the black circles the points r^ and Vj which are integrated in eq. (1^), the dashed line the function h{rij) 
and the oriented lines are related to the single particle wave functions. 

The density calculated in our model is correctly normalized. This can be seen considering that, since the uncorrelated 
term, corresponding to the diagram (a) of fig. ^, is already correctly normalized, the contribution to the normalization 
of the other terms should be zero. In the diagrammatic picture of fig. |l| the integration on ri transform in black the 
white dot. Since in the three point diagrams (d) and (e) of the figure the white circles are not reached by the 
correlation line hlra), the integration implies that the single particle wave functions connected to these points should 
be the same. Therefore the contribution of the diagram (d) cancel that of the diagram (6) as well as the contribution 
of the diagram (e) that of the diagram (c). 

We checked the validity of our model by comparing the density distributions we have obtained with those produced 
by a FHNC calculation. The comparison has been done using correlations and single particle wave functions given by 
the minimization of the energy functional generated by the semi-realistic Afnan and Tang S3 interaction H . 

In fig. g the correlations functions for the ^^C and '^^Ca nuclei are presented. The dashed lines show the results of 
a gaussian parameterization (two free parameters), while the full lines have been obtained with a procedure, we call 
it Euler minimization, consisting in minimizing the energy functional calculated up to the second order. In this last 
calculation there is only one free parameter: the healing distance. 

The comparison between the results of the FHNC calculation and those of our model is done in fig. 0. The one- 
body densities are compared in the upper panels for calculations done with both kinds of correlations. The agreement 
between the result of the two calculations is very good. Since the one-body densities are the diagonal part of the 
density matrix, in order to test our model also on the off diagonal part of this matrix, the comparison with the full 
calculation has been done also for the momentum distributions. These results are shown in the lower panels of fig. 0, 
and also in this case the agreement between the two calculations is rather satisfactory. 
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FIG. 1. Diagrams representing the terms considered in our model for the calculation of the density distribution. 
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FIG. 2. Correlation functions obtained for ^^C and ''**Ca using the semi-realistic S3 nucleon-nucleon interaction. Dashed 
lines are gaussian correlations and full lines are the Euler correlation functions. 
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FIG. 3. Comparison between densities (upper panels) and momentum distributions (lower panels) calculated with our model 
(dashed lines) and with a FHNC procedure (full lines). In all the panels the upper curves show the calculations done with 
Euler correlations, and the lower curves those done with gaussian correlations. 



III. THE MODEL: EXCITED STATES 



The description of nuclear responses in the quasi-elastic region is based upon the approach developed by Fantoni 
and Pandharipande |p| and applied to nuclear matter. The basic ansatz consists in assuming that the nuclear excited 
states can be described as a product of a Slater determinant and a many-body correlation function which is the same 
correlation function describing the ground state: 



l*/)=^l*/) 



(8) 



In eq. (g|), |$/) indicates a Slater determinant which differs from |$o) by the fact that a certain number of hole single 
particle states have been substituted with particle wave functions. 

The many-body nuclear response at energy lo and momentum transfer q to a an external field (5(q) can be written 
as: 
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where the sum is understood to run on all the nuclear excited states. The evaluation of cq. (|9|) requires the calculation 
of the full cluster expansion M . 

Like in the case of the density distribution, instead of performing the full expansion, after the elimination of the 
unlinked diagrams, we consider only the terms containing a single correlation function h{rij). 

In the case the nuclear final state is characterized by one particle above the Fermi surface, the matrix element to 
be calculated are of the kind: 



^ph(q) = 



{%h\F+Q+iq)F\^o) 



(10) 



where we have indicated with \^ph) the Slater determinant where the hole single particle wave function h has been 
substituted by the particle wave function p. 

If the Q operator is the charge operator, eq. (HG) represents, in the limit for q — > and p ^^ h, the normalization 
condition of the density. This fact allows us to identify the diagrams which we should consider to have a properly 
normalized many-body wave function. The diagrams we consider are shown in fig. m One can reconstruct the 
diagrams of fig. ^ by closing the particle and hole lines. 
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FIG. 4. Diagrams considered in the calculation of the one-particle one-hole response. The symbols have the same meaning 
of those in fig hi The black squares represent the point where the external field is acting on. 
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FIG. 5. In (a) the nuclear matter proton responses calculated with our model (dashed lines) are compared with those 
obtained with a full FHNC calculation (full lines). The lower panel show the difference between the responses calculated in the 
two approaches. 

We have compared the nuclear matter charge responses calculated with our model |Ig] with those obtained with 
the FHNC procedure. The two calculations have been done with the gaussian correlation function of ^^C of fig. ||. 
The comparison between the results of the two calculations is shown in fig. |g, where the dashed lines (our model) are 
almost exactly overlapped to the full ones (FHNC). 

The excellent agreement between the two calculations has been obtained because we have included in our model, 
in order to have a correctly normalized many-body wave function, both two and three body terms. The effect of the 
three point diagrams can be seen in the figs. [^ and 7 where longitudinal an transverse responses for infinite nuclear 
matter and for the ^^C nucleus are compared. We have used also for these calculations the ^^C gaussian correlation of 
fig. ^ The finite system responses have been evaluated with the single particle wave functions producing the density 
shown in fig. pi 

In fig. 0, the left panels show the nuclear matter longitudinal responses for three different values of the momentum 
transfer. The full lines are the uncorrelated responses, i.e. the Fermi gas responses. The dotted lines have been 
obtained by evaluating only the two-points diagrams, while the results of the calculations where both two and three 
points diagrams have been taken into account are shown by the dashed lines. The effect of the inclusion of the three 
point diagrams has opposite sign with respect to that of the two points diagrams. This fact is not surprising since, 
in the normalization of the charge, the contribution of the the three points diagrams cancels exactly that of the two 
points diagrams. The finite size effect do not modify this conclusion as one can observe in the three right panels of 
the figure. 

The calculation of the transverse response has been done by considering only the one-body magnetization current. 
The left panels of fig. 7 show the nuclear matter responses. Also in this case the effect of the three-points diagrams 
has opposite sign with respect to that of the two-points diagrams. We observe that the final response is extremely 



close to the uncorrelated response. The different behavior with respect to the case of the longitudinal response is 
produced by the fact that, for the nuclear matter transverse response, the diagrams (d) and (l) of fig. |4 do not 
contribute. This is strictly true only for the infinite system. In finite systems, only the diagram d is zero. The right 
panels of fig. 7 show that this is the main responsible for the different behavior of the correlations in the longitudinal 
and the transverse response. 

The model has been extended also to the case when the final state has two particles in the continuum |1^] , by using 
the same procedure described above for the Ip-lh excitation. In fig. g we show the longitudinal 2p-2h responses 
for three different values of the momentum transfer . The dashed lines have been obtained considering only the two 
points diagrams, while the full lines show the results obtained considering both two and three points diagrams. Also 
in this case the contribution of the three-points diagrams is relevant, close to the 40% for some value of the excitation 
energy. 

It is necessary to remark that the absolute value of the response is three order of magnitude smaller than the Ip-lh 
response. 
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FIG. 6. Longitudinal responses at various values of the momentum transfer calculated for Ip-lh excitation. The full lines 
represent the uncorrelated response the dotted lines the responses obtained by adding the two-points diagrams, and the dashed 
lined those obtained including also the three-points diagrams. The curves shown in the three left panels have been obtained 
for an infinite system while those shown in the right panels have been obtained for ^^C using the wave functions giving the 
densities of fig. H. 
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Figure 7. Same as in fig. g for tlie transverse responses. 



IV. SUMMARY AND CONCLUSIONS 

We have developed a model to describe nuclear responses taking into account short-range correlations. We have 
shown that mean values of the density operator can be calculated with minimal error by using a cluster expansion 
truncated up to the first order terms in the function h[rij). This truncation should be properly done to maintain 
the correct normalization of the wave function. This implies that both two and three points diagrams should be 
considered. We have shown that in both Ip-lh and 2p-2h responses the quantitative contribution of the three-point 
diagrams is relevant and cannot be simulated by multiplying the responses by a constant renormalization factor. 

The calculations we have presented are limited to inclusive responses, but we plan to apply the model to treat 
exclusive processes. 
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